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We predict large cooperative effect involving two atom two photon vacuum Rabi oscillations in 
a high quality cavity. The two photon emission occurs as a result of simultaneous de-excitation of 
both atoms with two photon resonance condition uii + ui2 ~ u a + Ub, where uii, u)2 are the atomic 
transition frequencies and u a , t±>b are the frequencies of the emitted photons. The actual resonance 
condition depends on the vacuum Rabi couplings. The effect can be realized either with identical 
atoms in a bimodal cavity or with nonidentical atoms in a single mode cavity. 

PACS numbers: 42.50.-p, 42.50.Fx, 42.50.Hz, 42.50.Pq 



I. INTRODUCTION 



High quality cavities have led to the study of a new regime of radiation matter interaction viz the study of strongly 
interacting systems. Several newphenomena such as vacuum Rabi splittings 0,0,0,0], collapse and revival of Rabi 
oscillations 5], trapping states [a], and systems like micromasers 0,01 have been studied. More recent applications 
on high quality cavities are in the context of quantum computation |9(. Most of these works concern the interaction 
of the individual atoms. Earlier cooperative effects like optical bistability involving a large number of atoms have 
been investigated |10| . A large part of these studies concerns the situations where the atomic transition frequency is 
almost equal to the cavity frequency. In this paper we report an unusual cooperative effect involving two atoms in 
a nonresonant cavity. This cooperative effect arises from the simultaneous de-excitation of two atoms such that the 
sum of the energies of emitted photons is equal to the sum of the excitation energies of the atoms. We demonstrate 
that in a high quality cavity the two atom two photon resonant effect could be large thus opening up the possibility of 
a variety of nonlinear i.e. multi-photon cooperative phenomena in nonresonant cavities. For this purpose the recent 
development on the trapping of atom inside the cavity should be especially useful. We bring out the origin of 
such large two atom two photon Rabi oscillations. 

We start by noting that in a two photon emission process the two photon resonance between the excited state |e) 
and the ground state \g) would occur at a frequency given by u> eg = 2u>, where w eg is the atomic transition frequency 
and lo is the frequency of the photons emitted. The process proceeds via intermediate states \i), which are away 
from a single photon resonance. Now consider an inter-atomic process involving two atoms with distinct transition 
frequencies lo\ and u> 2 such that lo\— lo and ll>2 — to are large so that individual emissions are not important. However, 
as shown in Fig^a), one can consider a two photon emission process such that lj\ + ll>2 — 2lo. Clearly this would be 
a cooperative process as it involves two atoms. Besides it should also be important as it is a resonant process. Let 
us then examine the transition probability for such a two photon emission. Let H + be the interaction responsible for 
the emission of a photon defined by the interaction Hamiltonian which is written in the form 



Hj 



(1) 



Then the second order perturbation theory leads to the following expression for the rate of two photon emission 



2n 



(ffi,S2|#+|si,e 2 )(5i,e2|-ff + |ei,e 2 ) , {g 1 ,g2\H + \e 1 ,g 2 ){ei,g2\H + \e 1 ,e 2 ) 



H(UJ2 — cj) 



5(u>l + UJ2 — 2u>). 



(2) 



Note that surprisingly R c = 0, as the two photon matrix element vanishes when oji + = 2lo as there are two paths 
for two photon emission which interfere destructively. It has been argued that a nonzero two photon emission can 
result if we include inter-atomic interactions [12j,ll3| which, however, are important only if the inter-atomic separation 
is less than a wavelength. A remarkable demonstration of such two photon cooperative effects is given in a recent 
work |lH using the methods of single molecule spectroscopy. Similar results apply to the case of two photon emission 
by identical atoms (Fig. 1(b)) if the photons of frequencies tu a and u>b are emitted 



LJ a +UJ b = 2w 



(3) 



In this paper we examine such two photon emission processes in a cavity. It is advantageous to use a cavity for 
the study of such a fundamental process as one would not be constrained by the requirement of small inter-atomic 
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FIG. 1: Two ways for two atom two photon emission, (a) corresponding to two possible intermediate states |ei, g-z) and |gi,e2) 
in the system of nonidentical atoms interacting with a single mode vacuum, (b) in the system of identical atoms interacting 
with two modes of the vacuum. 



separation. We demonstrate how high quality cavities can lead to a large two photon Rabi oscillation involving two 
atoms. Note that vacuum Rabi oscillations in the context of a single atom interacting strongly with vacuum inside a 
single mode resonant cavity have been studied extensively P, 0, S H • We also note that the two photon micromaser 
in a single mode cavity has been realized Q. In this paper we consider two different cases of two photon vacuum 
Rabi oscillations, (?) two identical atoms interacting with vacuum in a two mode cavity, (ii) two nonidentical atoms 
in a single mode cavity. 

The paper is organized as follow. In Sec. II we consider the case of two identical atoms interacting with two modes 
of a cavity and discuss two photon vacuum Rabi oscillations when the photons are emitted in different modes under 
resonance condition. In Sec. Ill wc consider the case of two nonidentical atoms interacting with a single mode of the 
cavity. We present both approximate and analytical results. In Sec. IV, we confirm that the two photon vacuum Rabi 
oscillations survive in the limit of small damping in a high quality cavity. Finally in Sec.V, we present conclusions 
and future outlook. 



II. TWO IDENTICAL ATOMS INTERACTING WITH VACUUM IN A BIMODAL CAVITY 

We consider two identical two level atoms, with transition frequency ujq, interacting with two modes of the vacuum 
having frequencies u> a and u>b in a cavity as shown in Fig^b). The Hamiltonian for the system is 

H = fno a a) a + fkdbwb + 

h [yd e *X e *l " \9i){9i\) + h}(ftl(.9ia + 32&) + \g t )(e t \(gi^ +92^)] , (4) 

i=l,2 

where a and (b and tf) are annihilation and creation operators for first(second) mode of the cavity, g\ and g 2 are 
the coupling constants. In a frame rotating with frequency loq, the Hamiltonian Q becomes 

H = -h^a-h5b%+^h[\e i ){g l \{gia + g 2 b) + \gi){e i \{g 1 a) +g 2 b^)\, (5) 

i=l,2 

A = U) - LU a , 5 = LU - LU b . 

We consider the special case of two photon emission i. e. the case when the initial state of the atom-cavity system is 

|V(0)) = | ei ,e 2 , 0,0). (6) 
Considering all possible states of the system in evolution, the state of the system at time t can be written as 

|V(t)> = ci(i)|ei, e 2 , 0, 0) + {\e u g 2 ) + \gx,e 2 )) {ca(t)|l, 0) + c 3 (t)|0, 1)} 

+\g 1 ,g 2 ){c i (t)\l, 1) + c 5 (t)\2, 0) + c 6 (£)|0, 2)}. (7) 
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Different terms in the wave function J7J) correspond to no photon emission, one photon emission and two photon 
emission. The photon emission can take place in either mode. A very interesting aspect of the state Q) is its 
entangled nature. This provides a method of producing entangled states, say, entanglement of two cavity modes pL5j . 
The time dependent amplitudes Ci(t) are determined by 

c'i = — «<7iv2c2 — ig 2 V2c3 

c 2 = iAc2 — igiV2c\ — ig 2 V2c4 — 2igic$ 

c 3 = iSc 3 - ig 2 \f2c x - ig 1 V2c4, - 2ig 2 c e 

c'4 = i(A + <5)c4 — ig 2 V2c 2 — ig\\[2c-i 

c 5 = 2iAc 5 - 2ig 1 c 2 

c' 6 = 2iSc 6 - 2ig 2 c 3 . (8) 

The complete solution of Eq.© has six eigenvalues corresponding to those there will be fifteen peak spectrum. In 
order to understand the nature of the two atom two photon resonance we present numerical as well as approximate 
analysis which can capture the physics of the cooperative process. We consider the case when detunings to the cavity 
field are much larger than the couplings i.e. |A| , \5\ >> g\ , g 2 but |A + <5| is small, the condition for two photon 
resonance is A + 5 = 0. In such a case cooperative two photon process should dominate and single photon processes 
would be insignificant. The results of numerical integration of Eq.JSJ are plotted in Fig|2 In the case when g\ 7^ g 2 




FIG. 2: (Color online) Two atom two photon emission probability, | C4 (t ) | 2 in a system of identical atoms interacting with 
vacuum in a two mode cavity, for (/2/gi = 1.5 and A/gi — —5.0. 

a novel resonance is achieved. The probability of two photon emission at resonance is quite high. The resonance is 
shifted from the position A + S = 0. This shift is due to the strong coupling to the vacuum field in the cavity. For 
gij g\ = 1-5 and A = — 5gi maximum two photon emission probability is approximately 0.9 and the interaction time 
required for achieving maximum probability is given by g\t « 67r. 

Having established numerically that the two photon resonance can be large in cavities, we present approximate 
analysis to demonstrate it. Under the above mentioned conditions for two photon resonance we can eliminate fast 
oscillating variables c 2l c 3 , c 5 , c 6 and effectively reduce the dynamics in terms of slowly oscillating variables C\ and C4. 
A simple treatment where one sets c 2 = C3 = C5 = cq = does not yield the physics of the two atom two photon 
emission. We thus relegate the procedure for eliminating fast variables to the appendix. The reduced form of the 
Eq.© is written as 

/ 2gjA 2gU \ / A S \ 

Cl = 1 Ia^M + s^l) Cl + 2 ^ { a^! + WTM) C4 
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/A 8 \ / c 2 5 f<5 2g|A , 



The solution of Eq. © gives 



M*)l 2 = i^r^ 8111 2 ■ (10) 

With G = 2 ^ ( + ^| ) ' ° = A + ' + 2(5? - ^ ( A^ - ^|) (ID 

Note that in the limit g\ = 92 and A + 8 = 0, the probability amplitude C4 for two photon emission tends to zero, as 
both Q and the numerator in Eg .1 (10)1 become proportional to (A + 8). Thus when couplings to the modes are same 
two photon emission probability has no resonance. In this case the transitions from |ei,e2,0,0) to \gi, g2,l,l) via 
states 4=(|ei,<72) + |<7i, e 2))|l, 0) and -^(|gi,e 2 ) + |ei , 52) ) |0, 1) interfere destructively. We further note that to order 
g\g\ the two photon resonance does not occur 

I /.M2 Iggigj • 2 St . 2 At 

l C4 (*)l = j 2 A 2 sm Y sm — ■ ( 12 ) 

The usual second order perturbation theory cannot lead to inter-atomic two photon resonance. One has to consider 
higher order terms in g\ and 52- However then the excitation itself would be negligible. Therefore one needs high 
quality cavities. The probability of cooperative emission of two photons in different modes is a periodic function of 




FIG. 3: The maximum value of the two atom two photon emission probability, |c4 (t) | 2 in the system of two identical atoms 
interacting with vacuum in a two mode cavity, is plotted with respect to (a) detuning 8 and (b) time, for gijgx = 1.5 and 
A = — lOgi. The solid line is corresponding to approximate result and the dotted line (...) corresponding to exact numerical 
result. 

time. In Fig|3{a), we plot the maximum value of | C4 (t) | 2 as a function of 8 and in Fig|3Jb) as a function of time t, for 
fixed values of g\, gi and A. At two photon resonance the probability corresponding to two photon emission in one of 
the two modes is much smaller than the probability of two photon emission in different modes. From Eqs. (|10|l and 
it is clear that the two photon resonance occurs at A + 8 + 4 (,g 2 /A + g\/8~) ~ 0. Thus the interaction with the 
cavity modifies the condition of two photon resonance. This is seen quite clearly in Fig|^a). We note the connection 
of the resonance frequency f2 to the one photon Stark shifts. It is well known that the shift in the frequency of a two 
level atom in the presence of a field with n photons is given by 2g 2 (n+ 1)/A which is equal to 4g 2 /A for n = 1. Thus 
the change 4 (g 2 /A + g^/o) is equal to the frequency shift of both the atoms due to the presence of a single photon. 
We have checked using the full solution of the Schrodinger equation that the result is quite good. However it 
should be borne in mind that the exact result is not periodic and exhibits rapid variations though the envelop agrees 
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with the result llOj l. The above mentioned approximate results are valid for larger values of detunings but for larger 
values of detutings a large interaction time is required to reach the maximum of two atom two photon transition 
probability. This should be possible with the recently developed method of trapping atoms in a cavity 0] . The other 
possibility is to work under the conditions of the Fig[3 



III. TWO PHOTON EMISSION BY TWO NONIDENTICAL ATOMS IN A SINGLE MODE CAVITY 

In this section we analyze a system of two nonidentical atoms interacting with a single mode vacuum field in a 
cavity (Fig|T[a)), Consider two nonidentical two level atoms having their excited states |ei), |e 2 ) and their ground 
states \gi), \g 2 ) interacting with a single mode cavity-field of frequency u. The Hamiltonian of this system is 



^(| ei )( ei | - \ 9l )( gi \) + -± (|e 2 )(e 2 | - \g 2 ){g 2 \) + coa^a 



H = h 

+h 9l (|ei)((?i|a + ct\gi)(ei\) + hg 2 (\e 2 ){g 2 \a + a t |<? 2 )(e 2 |) , (13) 



where uji(lu 2 ) is transition frequency for first (second) atom, a and are annihilation and creation operators for 
the field, and <?i(<? 2 ) is the coupling constant to the cavity mode with first(second) atom . In a rotating frame the 
Hamiltonian H can be written as 

H = -hA\g 1 )(g 1 \-h8\g 2 )(g 2 \ + hg 1 (\e 1 )(g 1 \a + a^\g 1 )(e 1 \) +hg 2 (\e 2 )(g 2 \a + a^\g 2 )(e 2 \) , (14) 
A = u>i — ui, S — uj 2 — lo. (15) 

Let us consider an initial state IV'(O)) — |ei,e 2 ,0) with both atoms in the excited state and cavity in the vacuum 
state. The state of the system at time t can be written as 

|V(t)) = c 1 (t)|e 1 ,e 2 ,0)+c 2 (t)|e 1 , 52 ,l) + c 3 (t)|gi,e 2) l)+C4(t)|pi,52,2), (16) 

where the expansion coefficients c's satisfy 

c'i = -ig 2 c 2 - igic 3 , 

c' 2 = iSc 2 — igiV^a — ig 2 ci, 

c'3 = iAc 3 — ig 2 V2c4 — igid, 

64 = i (A + S) c 4 — igiV2c 2 — ig 2 V~2c 3 . (17) 

The two photon resonance condition for this system would be A + S = 0. For couplings g\ , g 2 , much smaller than 
|A|, I S I, the solution of Eq.^TTJ) gives 

4g lg2 V2 . St . At 

C4(t) = — sin — sin — + higher order terms. (18) 

oA 2 2 

The first term in Eq. l|18|) represents independent emission by each atom. Clearly, to lowest order in gig 2 no two 
photon resonance occurs. Such a resonance can come from the terms of the higher order. Assuming that |A| and |<5| 
are large but |A + S\ is small, we eliminate fast oscillating variables c 2 and C3 in a way similar to the previous case 
and the Ea. (|17fl . in terms of slowly oscillating variables reduces, to 

(9x d\ ■ /7i( A 6 \ 

C 4 = i p lfls ,^(_^ + J: ^) Cl+i (A + 5 + M + M) C4 . (19) 
We find the approximate result for the two photon emission probability 



4G' 2 . 2 V4G' 2 + fl'H 

| C4 WI 2 = ^T^ W 2 < (20) 
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For large |A| and \S\ the Eq.JSUJ) shows two photon resonance at A + S + 3(gf/A + g^/S) ~ 0. Further such two atom 
two photon resonance appears for gi ^ g 2 , which disappears when gi — g 2 . In the latter case the antisymmetric state 
(|<7i,e2,l) — |ei,g2,l)) /V2 is decoupled from |ei,e2,0) and \gi,g 2 ,2). We present numerical results in Fig0] . The 
graph shows two photon resonance for gi ^ g 2 . It is clear that the position of resonance is shifted from A + 8 = 0. 
This shift in the position of resonance is due to larger values of g\ and 172, and depends on the ratio g 2 /gi- There is 
a large enhancement in the probability of two photon resonant emission in a high quality cavity. It is expected that 
such effects can be studied by placing the system used by Hettich et aZ0| m a cavity. 




FIG. 4: (Color online) Two atom two photon emission probability, |c4 (i) | 2 in a system of nonidentical atoms interacting with 
vacuum in a single mode cavity, for A/gi — —5.0 and gijgx = 2.0. 



IV. EFFECTS OF CAVITY DAMPING 

Before concluding we examine the effect of cavity decay on two atom two photon vacuum Rabi oscillations. We 
do a calculation based on master equation. Let 2n a and 2k;, be the rate of loss of photons from the first mode and 
the second mode respectively. The density matrix of the system of two atoms interacting with two mode field in the 
cavity will evolve according to the master equation 

p = —z[H, p] - n a (a'ap - 2apa) + pa) a) - K b (tfbp - 2bptf + ptfb) . (22) 

The density matrix for this system can be expressed in terms of all the states which are generated by the combined 
effect of H and dissipation. Because of the cavity decay, many more states are involved in the dynamics. For example 
for identical atoms interacting in a bimodal cavity, the relevant states are |ei,e2,0, 0), \g\, 0, 0), \g\, e 2 , 1, 0), 
lsi,e 2 ,0,l), |ei, 32, 0,0), \ei,g 2 , 1,0), |ei,g 2 ,0,l), |ffi,32,0,0), \g 1: g 2 , 0,1), \gi,g 2 ,l,0), \g!,g 2 , 0,2), \g 1: g 2 , 1,1), and 
\gi, g 2 , 2, 0). For this system density matrix is expressed as 

1 i'+j' i+j i'+j'-k' i+j-k 

p= E EE E p(*',j',k',i f ,t,j,k,i)\i',f,k',i')(2,j,k,i\ . (23) 

i',j',i,j=0 fc'=0 fc=0 l'=0 1=0 

Here i, i' (j,j f ) represent states of the first (second) atom with the convention |0) corresponding to excited state and 
|1) corresponding to ground state, the indices k, k' (I, V) represent the number of photons in the first (second) mode. 
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Thus the dissipation requires considerable numerical work. Results for two identical atoms in a bimodal cavity are 
shown in Fig[SJ We show results for optical cavities with g/n 30 in FigJSfb) and for currently realizable cavities 
(g/n = 10) in FigE^c). The two atom two photon vacuum Rabi oscillations survive in the limit of small damping 
g/n ~ 30 but for larger damping {g/n — 10) die fast. Similar results are found for two nonidentical atoms in a single 
mode cavity. 




FIG. 5: Periodic behavior of two atom two photon emission probability | C4 (i ) | 2 , for identical atoms interacting with vacuum in 
a bimodal cavity, for S = 3.5gi, A = — 5<?i, Q2 = L5gi and cavity damping constants (a) K a = Kb = 0.00, (6) K a = Kb = 0.03(?i, 
(c) Ka = Kb = O.lffl. 



V. CONCLUSIONS 



We have reported large two atom two photon vacuum Rabi oscillations in two systems, one having two identical 
atoms in a two-mode cavity and another having two nonidentical atoms in a single-mode cavity. We have shown that 
for asymmetric couplings {g\ 7^ 52), the probability of two photon emission is quite large but for symmetric couplings 
(.9i = 92)1 the two photon emission probability is very small. Further, we have shown that the condition of two 
photon resonance in the case of strong atom- field interaction is modified from its free-space form (A + S = 0). These 
two photon transitions involving two atoms can be used for generating and detecting different types of entanglement 
between two field modes and two atoms 

GSA thanks V. Sandoghdar and G. Rempe for discussions on the two atom two photon resonance. 
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APPENDIX A 



Our procedure for eliminating fast oscillating variables is extended form of the procedure discussed in Ref. [TtI ] . The 
Hamiltonian JSJ can be written as 



H = Hq + eV, 
where 

H = -Aa^a-Sb^b, eV = ^ h [\e t ) ( gi \( gi a + g 2 b) + \g l )(e t \(g 1 a t + g 2 b i )] 



(Al) 



The eigenstates and corresponding eigenvalues of Hq are 

1) = |ei,e 2 ,0,0) 

2) =2- 1 /2 ( | ei;ff2) + | fflje2)) | li o) 

3) =2- 1 / 2 (|e 1 , 52 ) + | 5l ,e 2 ))|0,l) 

4) = 1^,52, 1,1) 

5) = |.gi,.92,2,0) 

6) = | 5 i,52,0,2) 



Ei = 0, 
E 2 = -A, 
E 3 = -S, 
E A = -(A + S), 
E 5 = -2A, 
E 6 = -26. 



The resolvent for Hq is the function 



G {z) 



1 



z-Hq' 

where z is complex. If Pj is projection operator for the eigenstates of H 

Pi = \i){i\; i = l,2...6. 

The resolvent Go can be expressed as 

Pi 



z-E, 



The resolvent for the full Hamiltonian H is 

G{z) 



1 



z-Ha-eV' 
1 



1 + eV 



1 



z-h r 



z-H 
= G (l + eVG). 

From Ea . (|A5|I the resolvent for the full Hamiltonian H can be expressed in the power series of e as 



(A2) 



(A3) 



(A4) 



G = ^e n Go(l/G )". 



(A5) 



(A6) 



For small values of e, G(z) has singularities in the complex z-plane in the neighborhood of poles of function Go i.e. 
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FIG. 6: The contour in complex plane shielding two eigenvalues -Ei and E& and leaving others outside. 



eigenvalues of Hq. Further eigenvalues E\ and E 4 are very close to each other under the condition A + S « and 
other eigenvalues are largely separated. We consider a contour, T in the z-plane that encloses eigenvalues E\ and E 4 
only and leaves others outside as shown in the FigEl We define a new projection operator Pp as 



P r = Pi+Pi 
1 



2iir j T 



G{z)dz. (A7) 



Here Pi and P4 are the projection operators for eigenstates of full Hamiltonian H corresponding to the eigenvalues 
inside the contour. The effective Hamiltonian will have the form 

H eff = (P 1 + P i )HP v (P 1 + P 4 ). (A8) 

From the definition of the resolvent we have 

(z — H)G = G{z — H) = 1. 

HP r = — I zG(z)dz. (A9) 
1m J r 

Substituting value of G(z) from Eg. (|A6|) in Eg. (|A9|> and interchanging summation to the integration we have 

HP r = V ^- I zGo(VG ) n dz. (A10) 

The effective Hamiltonian can be expressed as 

00 

H eff = E 1 P 1 +E 4 P 4 + ^2e n A ( - n '>; 

n=l 

00 

A (n) = (Pi + p 4) J- I zG {VG a ) n dz{P l + P 4 ). (All) 

' — ! Z17T r 
n=l JL 

Inside the contour T, Go has singularities at E\ and E 4 only so the integral in the Eci. (|All|l is nothing but the sum 
of the residues at z = E\ and z = E 4 . Further as in our case ePiVPi, eP 4 VP 4 and eP\VP 4 equal to zero, there is no 
first order and third order terms. The second order term is 

A™ = PxVQxVPx + P 4 VQ 4 VP 4 + PxVQiVPi + P 4 VQ 4 VP X ; (A12) 



The forth order term is 



2iwJ T \z-Ex z~E A ) gf^z-Ei \z-Ex z - E 4 ^z-Ejj 

y _^L_y ( + ^L_) dz (A13) 
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For simplification we use the condition for resonance A + S — 0, i.e. E\ — E4. Thus the forth order term is 



A (4) 



1 



Pi 



Pi 



z- E 1 



Ex 



V 



P 



^— ' z 



J — u V - 

£7,- ^ z 



E, 



-V 



Pi 



Pi 



Ei 



Ei 



Integrating Ea. (jA14j) we have the forth order term 

= (Pi + VQiVQiVQiViPi + Pi). 



dz. (A14) 



(A15) 



Using the values of Ei, E2, E3, Ei, E$, Eq and V the effective Hamiltonian expressed in basis |ei,e2,0, 0) and 
(51,52, 1, 1) is 



H eff ~ 



— + — - 

2ffl32 _ 23132 
A 5 



, 4 9l 4g 2 

1 4g 1 g 2 , igig 2 
"T A 3 + <5 3 



2flig 2 2gig 2 , 4g^g 2 , 4gigf 
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(A16) 



With some algebraic manipulation and considering g\ and 52 up to forth order effectively the Hamiltonian (0 reduces 
to 



Heff = 
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2g^ 



(A17) 



It should be noted here as two atom two photon resonance appears at large interaction time in dispersive limit, the 
terms in the effective Hamiltonian up to forth order are important to predict correct evolution. Using the effective 
Hamiltonian HA17f) the Eq. © reduces to Eq. © . 



